Quantum stochastic description of collisions in a canonical Bose gas 
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We derive a stochastic process that describes the kinetics of a one-dimensional Bose gas in a 
regime where three body collisions are important. In this situation the system becomes non inte- 
grable offering the possibility to investigate dissipative phenomena more simply compared to higher 
dimensional gases. Unlike the quantum Boltzmann equation describing the average momentum 
distribution, the stochastic approach allows a description of higher-order correlation functions in a 
canonical ensemble. As will be shown, this ensemble differs drastically from the grand canonical 
one. We illustrate the use of this method by determining the time evolution of the momentum mode 
particle number distribution and the static structure factor during the evaporative cooling process. 

PACS numbers: 03.75.-b,03.75.Lm,05.45.Yv 
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The stochastic kinetic description of a quantum gas 
has been the subject of considerable activity [1-3]. The 
main idea is to describe the state in terms of a set of 
stochastic classical fields tp^'(r,i) (i = 1, ...,/) associ- 
ated to the bosonic quantum field ip{r, t) whose averages 
give access to the n-point correlation functions, possibly 
in the CN (canonical) ensemble. The determination of 
these quantities has been shown to be of relevance in re- 
cent interference experiments with cold atoms Q where 
mode probability functions have triggered considerable 
interest 0, ||. As discussed in [l|, the main difficulty 
is to be able to derive an associated Fokker-Planck(FP)- 
like equation for the probability in the P representation 
which can be recasted into a stochastic Langevin equa- 
tion. Such a task cannot be achieved exactly since the 
resulting equation would be a non dissipative FP equa- 
tion [JH-i|. Only the use of the positive P representa- 
tion appears to be tractable @- Some approximations 
have to be made in a bid to obtain numerically solvable 
dissipative equations. For the general case of finite tem- 
perature, only a phenomenological approach is available 
in which the stochastic terms guarantee the approach to- 
wards thermodynamic equilibrium without the require- 
ment of an additional cutoff at large momentum [H . 

In this Letter, we present a microscopic derivation of 
a Langevin equation starting from the full many body 
Hamiltonian in which the stochasticity originates from 
collision processes. For simplicity, we study a one dimen- 
sional Bose gas subject to three-body interactions only 
(two-body interactions does not lead to effective ther- 
malization in ID). Such a system could be realized in 
cold atom experiments fioj and differs from the usual ID 
exactly solvable Hamiltonian by the addition to the two- 
body quartic term of a sixtic interaction term resulting 
from virtual interaction with transverse degrees of free- 
dom. Its study allows simpler access to the exploration of 
quantum collisions com par ed to the more involved higher 
dimensional systems (3. 1 111 Il2j , in particular to the study 
of their enhancement due to the Bose statistics. As an il- 
lustration, we consider the dynamics of evaporative cool- 
ing in the CN ensemble and show that the time evolution 



of quantities such as correlations or mode probability dis- 
tribution differ strongly from the one obtained from the 
Boltzmann equation valid only in the grand canonical 
(GN) ensemble. 

Defining the atom mass m, the scattering length a s , 
its ID density n and the transverse trap frequency uj r , we 
focus on the thermalization regime described in [10( char- 
acterized by a high transverse confinement Tiuj r 3> fc^T, 
a weak correlation parameter 7 = 2a s mui r /Tin < 1 in 
order to violate integrability so that the three-body col- 
lision rate overcomes the two-body collision rate respon- 
sible for transverse mode excitations. 

Restricting ourselves to a uniform gas and defining 
the spatial Fourier components of the field tp^ (r, t) — 
exp(ikr)a k {t)/y/L, we shall show how to derive the 
following Langevin-type equation: 
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where uj' p is the kinetic and mean field energy of the atom, 
r* n (f) and T° ut (t) correspond to the ingoing and outgo- 
ing collision terms respectively. These are averaged func- 
tionals of the a p (t) and are determined from Eq. (|17ll8p . 
The noise dr]^ follows a Gaussian distribution with the 
only non trivial average (dr) p dr]*i) = S p ydt. In contrast 
to @, 01 , this equation includes both high and low energy 
modes in the stochastic process. The average is defined 
as (A) = ^2l =1 / 1 for any set of realization AW. It 
is done in the grand canonical (GC) ensemble and cor- 
responds to the Monte-Carlo approximation of the inte- 
gral over coherent state labeled by the a p 's @. In other 
words, the expectation value of any observable functional 
of the creation annihilation operator a p and a p , is in the 
limit of a large set /, identical to an average over the 
stochastic variables a* and a p . 

We note the connection with the quantum Boltzmann 
equation. Multiplying ([1]) by its complex conjugate and 
taking the average over the ensemble, the stochastic 
equation is connected to the quantum Boltzmann equa- 
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(Kl 2 ): 

2Y k n {l + n k )-2T° k ut n k (2) 



tion for the momentum distribution n k — (n k 

dn k {t) 
dt 



This equation becomes closed if the T k n and r£ ut are 
functionals of rl k . This is realized through the applica- 
tion of the Wick's decomposition theorem (which cor- 
responds to the stosszahlansatz): defining the product 
M n ({a k }) = a* ki . . . a* k a k ^ . . . a k / n , its average in the GC 
ensemble for a uniform gas is decomposed as: 



(M n ({a k })) = n kl . ..n k , 
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where P is the permutation ensemble. On the contrary, 
the stosszahlansatz is not necessary anymore in Eq.([T]) 
and correlations can be taken into account. 

Another advantage of the stochastic formalism is the 
possibility to relate the averages in the CN ensem- 
ble to the GC ones by means of the weight functions 



W N ({a k }) = e ~£>> 
atoms we obtain for any product 
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(W N ({a k })) 
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The weight function plays the role of a projection opera- 
tor restricting the particle number to N. Probability dis- 
tribution in the CN ensemble can also be determined for 
mode population. For example, defining the weight func- 
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tions WN,kjto 

W N , k=0 = e-M 2 \a Q 
for the zero momentum mode reads: 



• ^ u|Qt ' 2 (£^oKI 2 J /N\ and 
2N /Nl, the probability distribution 
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The determination of the stochastic equation is done as 
follows. We start from the ternary many body Hamilto- 
nian Qjjl: 



H = J2 uj > 
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where uj k = k 2 /2m is the atom kinetic energy and 
gs = 2\n{A/i)huj r a 2 is the ternary three body interac- 
tions. The suppression of the quartic interactions pro- 
portional to g 2 = huj r a s is valid in a region where the 
mean field interaction energy is negligible compared to 
the kinetic energy i.e. g 2 n -C fefiT so that the phonon 
excitation energy are mostly particle- like [Toj . Starting 
the density matrix pit), we obtain the reduced density 
matrix for the momentum mode p by taking the trace 
over all other modes: 



p p (t) = Tr\ p (p(t)) 



(7) 



The stochastic time evolution of the mode p is derived by 
considering the other modes as a bath for this mode. For 
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FIG. 1. Evolution of the momentum average distribution nk 
vs. the reduced time r. The corresponding entropy is shown 
below and indicate how the distribution is close to equilib- 
rium. Initial and final distributions are shown in blue and 
red respectively in the inset. 



this purpose, we decompose the density matrix in terms 
of an uncorrclatcd contribution and a correlated one: 



p(t) = Y[p P {t) + 6 P (t) 



(8) 



The Hamiltonian can be decomposed into three terms: 

H = H 1 +H 2 + H 3 (9) 
where n k = a k a k and 

H 1= uj p n p + 0(1/L) (10) 

H 2 = A p n p + Cla p + alC p + 0{l/L) (11) 

and where A pi C p and H3 are operators describing the 
other modes. We omit terms of the order 0(1/ L) <C n 
which are negligible in the thermodynamic limit. Using 
the formalism in [H to derive the master equation up to 
the second order in g%, we obtain: 

^ = -^[H[,p p }+r p n (t)(2alp p a p -[a p alp p } + ) 

- V; ut (t)([ala p , p p ] + - 2a pPp al) (12) 

where 

H[=H 1+ Tr\ p (JJ p k H 2 ) = J p n p (13) 
k^p 

defining oj' p = cu p — Slg^n 2 and n = ^2 k (n k } / L and where 
1 f°° 

T™(t)=Re- 2 - dt'(Cl(t')C p ) 
n Jo 

1 f°° 

= Re-^ / dt>Tr(Cl(t>)C p J[ p k (t)) (14) 

"> JO U~L„ 
1 f°° 

T out {t)=R dt' (C p {t')Cl) 

n Jo 



k^P 



(15) 
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Terms containing the operator A p do not contribute up 
to g 2 by symmetry between the in and the out terms. 
The only contribution comes from the operator C p which 
expressed in the interaction picture becomes: 

CpJS) = ~~j~2~ e P * ^ C P2 C P3 C 9i C 92 C 93(16) 

I2i Pi=I2i m 

The stochastic equation can be implemented numerically 
I 



directly using Eq. (|14ll5l) but in order to avoid the sum- 
mation over too many momentum variables we make the 
random phase approximation. For a uniform gas, we can 
neglect off-diagonal contributions in the pi and qi since 
they induce a phase factor that appears to be random 
and thus cancels in the summation process. Thus, only 
diagonal components remain. Carrying out the integral 
over t' , we obtain with the following expression: 



I™(t) = (3!) 2 tt (jgL) 5 E lP .,E, 9 ^(E^-E^)K^f3K + l)K 2 + l)K 3 + l)> (17) 

{Pi,qi} « « 

Note the Bose enhancement factor (terms n p + 1) that amplifies the collision process when the output modes are 
already populated. From this form we can deduce the FP equation associated to the master equation (fT2")l : 

= W " + (rr (t) " r >" (<)) (^ ap + + 2r " w ^4} Pp(<) (19) 



from which we deduce the Langevin equation Eq.(JTJ. 
The solution of this equation is a Gaussian distribution 
P p (t) — cxp(— \a p \ 2 /n p (t)) /n p (t) with the normalization 
/ d 2 a p P p (t)/n = 1 and where n p (t) has to fulfill Eq.(J2|). 
In equilibrium in the GC ensemble, through the applica- 
tion of Eq.© on Eq. (|17ll8[) . we recover the Bose-Einstein 
distribution rik = l/[exp((w p — /i)/fcsT) — 1] as the sta- 
tionary solution where the parameter [i defines the chem- 
ical potential. 

For an illustration of all these considerations, we ap- 
ply the stochastic method to the process of evaporative 
cooling. We take a gas of 87 Rb of n — 5.25 [im^ 1 con- 
fined in a box of size L = 20^m with u: r /2ir = 6kHz 
and a s = 5.3nm. We choose the dimensionlcss time 
t = lOH/t* where 1/t* = (3!) 3 (g 3 /hL) 2 m/nh. 

We start from an initial Bose gas at high tempera- 
ture Tj = 29nK, we then remove the hottest atom with 
momentum higher than k c = A^im" 1 and study the re- 
laxation process of this cut distribution towards an equi- 
librium one with a lower temperature Tj = 9nK at a 
time estimated to t = 0.14s (r = 200). The time evolu- 
tion of the GC distribution calculated from Boltzmann 
approach of Eq.® is shown in Fig.l together with the 
entropy S(t) = J2 p (™p + 1) l°g("p + 1) ~%> log(n P ) evolu- 
tion that allows to monitor the speed at which the equi- 
librium state is reached. Its production can be shown to 
be always positive and stops at equilibrium (l3j . 

The stochastic approach allows a more refined descrip- 
tion in the GC ensemble. For instance, the atom proba- 
bility distribution for the momentum mode is determined 
from the solution of Eq. p^)) and corresponds to a Poisson 




FIG. 2. Evolution of the CN and GC mean atom number in 
the zero mode and their fluctuations vs. the reduced time r. 
The negative correlation with the first excited mode is shown 
in the inset. 



distribution: 

1 P ' J 7T n p \ q (n p + l)" P + lV ' 

with the larger fluctuations {5 2 no)GC = no(na + 
1). Correlations between modes are non existent e.g. 
(SnoSrii) gc = where n\ is the mode for the first excited 
state k = 2ir/L. 

The situation changes drastically in a canonical de- 
scription as the fixed atom number restricts the possi- 
bility of fluctuations. Using the formulae (Q| and ([5]) 
over a sample of 10 6 variables for the stochastic process 
(HJ), if the average value remains practically unchanged, 
the fluctuations arc significantly reduced as shown as in 
Fig. 2. The process of evaporative cooling has the effect to 
increase these fluctuations together with the mean popu- 
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FIG. 3. Evolution of the CN atom probability distribution 
function in the zero mode vs. the reduced time r. For com- 
parison, the distribution for the mode k c is represented in the 
inset. 




FIG. 4. Evolution of the static structure factor S q in the CN 
ensemble. For comparison, the corresponding GC function 
Sg(t = 0) is represented in the inset. 



lation. The negative values of (SnoSm) c n is an evidence 



of correlations between modes in the CN ensemble. 

The canonical probability distribution resembles a 
Gaussian one. Fig. 3 shows how this distribution evolves 
with time and illustrates the benefit of using a stochastic 
approach for quantities that differ according ensembles. 
For the mode k c , the distribution is initially peaked and 
becomes monotonic during the time evolution. 

The static structure factor is defined in GC and CN 
ensembles respectively as fl4| : 

S q = (S Pq 8p_ q )/N Sf) = {Sp q 6 P - q ) N /N (21) 

and corresponds to the density fluctuations p q = 
Sfe c i- 9 c fe- These functions can be measured in Bragg 
spectroscopy experiments (l5| and reach unity for large 
momentum q. In FigUl we determine the time evolution 
of these fluctuations during the evaporative cooling pro- 
cess. We notice that the presence of correlation in the CN 
ensemble reduces significantly the static structure factor. 

In conclusion, we have reformulated kinetic theory in 
terms of a stochastic approach where the stochastic vari- 
ables describe the various momentum modes of the gas. 
Through a specific weight average, the CN description 
becomes accessible and shows drastic differences in com- 
parison with the GC one. The method has been used 
for a one-dimensional gas but is quite general and can 
be applied to more general cases of higher dimension, in- 
cluding higher order interaction terms or in presence of 
longitudinal trap confinement. 

PN would like to thank S. Heller, H. Hauptmann, S. 
Kronkc and W. Strunz for helpful discussions. 
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